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ABSTRACT 

An abstract approximation framework for the solution of operator algebraic Riccati 
equations is developed. The approach taken is based upon a formulation of the Riccati 
equation as an abstract nonlinear operator equation on the space of Hilbert-Schmidt opera- 
tors. Hilbert-Schmidt norm convergence of solutions to generic finite dimensional Galerkin 
approximations to the Riccati equation to the solution of the original infinite dimensional 
problem is argued. The application of the general theory is illustrated via an operator Ric- 
cati equation arising in the linear-quadratic design of an optimal feedback control law for a 
one dimensional heat/diffusion equation. Numerical results demonstrating the convergence 
of the associated Hilbert-Schmidt kernels are included. 
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ON HILBERT-SCHMIDT NORM CONVERGENCE OF 
GALERKIN APPROXIMATION FOR OPERATOR 
RICCATI EQUATIONS 

I.G. Rosen 


Abstract. An abstract approximation framework for the solution of operator algebraic 
Riccati equations is developed. The approach taken is based upon a formulation of the 
Riccati equation as an abstract nonlinear operator equation on the space of Hilbert- 
Schmidt operators. Hilbert-Schmidt norm convergence of solutions to generic finite 
dimensional Galerkin approximations to the Riccati equation to the solution of the 
original infinite dimensional problem is argued. The application of the general theory 
is illustrated via an operator Riccati equation arising in the linear-quadratic design 
of an optimal feedback control law for a one dimensional heat/diffusion equation. 
Numerical results demonstrating the convergence of the associated Hilbert-Schmidt 
kernels are included. 


1. Introduction 

In this paper we develop an abstract approximation theory for algebraic Riccati 
equations on spaces of Hilbert-Schmidt operators. Our approach is based upon 
Barbu’s [1] formulation of a class of Riccati equations as abstract nonlinear operator 
equations on a space of Hilbert-Schmidt operators. We argue that solutions to 
generic finite dimensional Galerkin approximations to the Riccati equation converge 
in Hilbert-Schmidt norm to the solution of the original infinite dimensional equation. 

Our effort here is closely related to results in one of our earlier papers [6] wherein 
we developed an approximation theory for operator Riccati differential equations 
using techniques similar to those which will be employed below. Our treatments 
here and in [6] differ from the standard approach to the analysis of operator Ric- 
cati equation approximation in that we consider the nonlinear operator equations 
directly in the space of Hilbert-Schmidt operators rather than integral equation 
equivalents and their limiting properties as the time horizon tends to infinity (see, 
for example, [3]). While we do in fact obtain a stronger convergence result than the 
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We seek a solution II e C to (2.3). We note that for n a positive integer, the 
operator B given by B($) = <$", $ eC (i.e.f(z) = z n ) can be shown to satisfy the 
conditions above. Indeed, boundedness follows from the esitmate < |$|^, 

while monotonicity can be established via 

n 

[$ n - ^ n ,$- > o, 

y= i 

for 4>, 4/ e C (see [9]). Finally, for 'k e C and A ^ 0 let be tbe orthonormal sot 

of eigenvectors of ^ with corresponding eigenvalues {o^}^. (The fact that $ 6 C 
implies oti > 0, i = 1,2,...). Then, if we let 7 ^ denote a nonnegative solution to 
the equation 7 ; + A 7 ™ — a; = 0 (the intermediate value theorem guarantees that 

00 

such a 7 ,■ exists with 0 < 7 ,- < a,) and set ^ 7 i(Vb v? e H, it follows 

1— 1 

that $ e C and $ + A$ n = \&. This establishes (2.2). When n = 2, (2.3) becomes 
the familiar quadratic Riccati equation. 

Define the operator A e £(V, V*) by 

+ $ e V. 

It can be argued (see [ 1 ]) that A is strongly V— elliptic - that is there exists a 
constant uj > 0 for which 

(2.4) [A*,*]n >u\M 2 v , $eV. 

If we define the subspace Dom(A) = {$ e V : A$ e Ti] then it follows (see [ 8 ]) 
that the operator A : Dom(A) C 7i — * 7{ is densely defined and m-accretive in 
'H. With the above definitions the problem of finding a solution to the operator 
algebraic Riccati equation (2.3) becomes one of finding a solution II e Dom{T) to 
the abstract nonlinear operator equation in H given by 

(2.5) F{Ti) = 0 

where 0 e C is given and T : Dom (F) C H — » Ti is the operator defined by 
:F($) = A§ + B($) for $ e Dom(F) =C D Dom(A). 

Using a standard fixed point argument on the closed convex subset C , Barbu [1] 
argues that the equation 

J\($a) = 0 

has a unique solution 4 >a € Dom(F) for each A > 0 where T\ : Dom(^F) C Ti —* Ti is 
the Yosida -like approximation to T given by T\ = A + A -1 {/ — (/ + AB) -1 }. Then 
using the boundedness and montonicity of B , Barbu argues further that the II a 
converge in 7i as A — ► 0 to an opeator II e Dom(^F) which is a solution to (2.5) (or, 
equivalently, (2.3)). The strong V-ellipticity of A (i.e. (2.4)), and the monotonicity 
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of B (i.e. (2.1)) are then used in the usual way to establish the uniqueness of II. 
Note that II e Dom (P) implies that II e C, i.e. that it is a nonnegative, self-adjoint 
opeator in "H, and that II e V with ,411 = A*II + IIA e Ti. 

3. Galerkin Approximation and Convergence Theory 

For each n = 1,2. ..let H n be a finite dimensional subspace of H with H n C V , 
n=l,2,... . Let P n : H — > H n be the orthogonal projection of H onto H n with 
respect to the (•, •) inner product on H. We assume 

(3.1) lim \\PnP ~ <p\\ = 0, ip e V. 

n— ►oo 

Note that assumption (3.1) implies that lim \P n p — p\ = 0, p e H, and that the 

n— *o o 

P n are uniformly bounded in the uniform operator topologies on C(H) and £(V). 

Lemma 3.1 The operators P n admit extensions, which we shall again call P„, to 
idempotent, uniformly bounded operators in £(V*). Moreover, lim || P n p — </?|| + = 

n — ►oo 

0 ,peV*, the V*-adjoint (i.e. the operator P* satisfying < P n p, >*=< ¥>■> P £ >* 
e V *) is given by P* = jPnl' 1 , and lim \\P*p - tp\\ m = 0,p e V*. 

n — hx) 

Proof For p e V* set P n p = p n where p n is the representer of the functional on 
H n which is the restriction of p to H n . That is (p n ,tp n ) = (P n y>, ip n ) = (c p,tp n ), 
e H n . It is clear that P n as given above is a well defined extension of the 
orthogonal projection of H onto H n and that it is idempotent. Moreover, since 
HnCV C V*, for p e V* we may consider P n p = p n e H n an element in V* 
via the duality pairing (P n p,xji) = (</?„,?/>), V* e V. Then for p e V*,ip e V, and 
<Pn = Pn<P we have (P n p,^) = (p n ,%l>) = (< pn,P n = (<p,P n ip)- Consequently for 
ip e V* we have 

ll-Pn^ll. = SU P l( P n^V’)| = Sup \(p,P n 1p)\ < ||y>||*||P B ||, 

ipeV iptV 

or 1 1 P n 1 1 * < ||P n ||. Thus assumption (3.1) implies that the P n are uniformly bounded 
in C(V*). (We note that alternatively the same extension of the projections P n to 
operators on V* could have been obtained by using the standard approach based 
upon the density of H in V* and the usual extension construction in terms of limits.) 
Now for p e H we have 

< P n p,il> >*= (P n p, 7 _1 V0 = (‘rSPnT -1 */’) = (<£>,7 -1 7P„7 _1 V>) 

=< v?,7P„7 -1 ^ >,, 

and consequently that the V*-adjoint of P n , P * , is given by P* = 'yP n j~ 1 . Finally, 
assumption (3.1) yields 

lim ||P*^-V?||*= lim ^l~ l P*p - 7“V|| 

n — >00 n— *oo " " 
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= lirn^ ||7 J 7P n 7 V“7 V|| 

= lim ||P„7“V “ 7~ V|| = 0 

n — >00 

for each <p e V* and the lemma is proved. 

Define the sequence of finite dimensional subspaces W n ,n = 1, 2, ... of H by 

W n = {$„P n :$„e£(P„)}. 

Clearly H n finite dimensional implies that all operators in H n are of finite rank and 
thus that H n is a subspace of both H and V. For each n = 1,2,... let C n C Ti n 
be the closed convex cone given by C n = {$ n Pn e H n • = ^n, > 0}. 

Note that C n C C,n = 1,2,... We define Galerkin approximations to the operator 
: Dom(F n ) C H n — ► H n , as follows: 

(3.2) P„($„P„) = {^$„Pn + P(^nPn)}|«„, for $„P„ e Dom(p n ) = C n . 

That is, for <fr„P„ e C„, .P n ($ n P n ) = *knPn € where ^nPn is that element 
in H n (guaranteed to exist and be unique by the Riesz Representation Theorem 
applied to the functional A$ n P n + P($ n P n ) e V* restricted to a functional on the 
finite dimensional Hilbert space H n ) which satisfies [A$ n P n + P($„P„), X n P n ]n — 
VnP mXnP h]hi XnPn^H n - 

It is of some value to note that the approximations to T given in (3.2) are the 
same ones that would be obtained via the standard approach which is based upon 
the replacement of the operators A and B in (2.3) by their respective Galerkin 
approximations on H n and H n . Indeed, for each n = 1,2,... define the operators 
A n e C(H n ) and B n : C n C H n — * H n by A n (p n = ip n , where for <p n e H n , if> n is that 
element in H n which satisfies («4y?„,Xn) = (V\i,Xn), X»» e H n , and B n (Q n P n ) = 
P n B(Q n P n )P n . From the fact that P„ is the orthogonal projection of H onto H n 
it follows that [$P n , = [$,'Ji r n ] 7 < for every $ e V* and \& n e H n . Then, for 

$„P„, ^„P„ e H n we have 

[Pn(^nPn), ^ nP n ]n = [A$ n P n + P($ n P n ), \I'„Pn]w 

= [{A* n P n + P(* n P„)}P n , <H n P n ]n 
= [A*$ n P„ + $«P„AP n + B($ n Pn)Pn, * „P«]w 

OO 

- E {(A*$ n P n e j fc,^ n P n e fe ) + ($ n P n AP n e Jfc ,'I'„P n e fc ) 

k = 1 

+(B($ n P n )P n e*, y n P n e k )} 

OO 

= E {(An$ n Pnek,'&nP n ek) +(A n P n e k , $* <b n P n e k ) 

k~l 

+(P„P($ n P„)P„e fc , ^„P„ejt)} 

OO 

— E ({^n^n + + B n ($ n Pn) } Pn^ki ^nPn^k)i 

k=l 
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or, 

Fn(*nPn) = {A* n $ n + $„A„ + B n ($ n P n )} P n . 

In the particular case when Z?($) = $ 2 , for example, the operators T n take the 
form Pni^nPn) = {A£$ n + $ n A n + $ 2 } P n . 

Set 0 n = P n QP n e C n and consider the problem of finding a solution II n e C n to 
the nonlinear operator equation 

(3.3) JF n (II n ) = ©n 

in 7i n . Arguments similar to those described in section 2 above yield that for each n 
= 1,2,..., the equaiton (3.3) admits a unique solution II n e C n . We shall argue that 
lim II !!„ — II||v = 0; that is, that the II n converge in the V-Hilbert-Schmidt norm 

n— ► oo 

to the solution II to the equation (2.5) (or equivalently (2.3)). In order to do this 
we shall require the following lemmas. 

Lemma 3.2 If {a,}^ is an absolutely summable sequence of real numbers, then 
there exist sequences and {c,}^, for which lim bi = 0, { c 2 - } , is absolutely 

t— *oo 

summable, and a; = i = 1,2, .... 

oo 

Proof Set a = ^ l a «l an d, for j = 0,1,2,... define the nonnegative integers kj by 
*=1 

k 0 = 0 and for j = 1,2,... let kj be the first index for which 

£>!>«-£• 

»=1 J 

Then setting bi = 1/j and c* = jai, for i — kj- 1 + 1, ...,kj, j = 1,2,..., we have 
b{Ci = a;, i = 1,2,..., lim b( = 0, and 

t— ►OO 

OO OO kj OO f kj kj - 1 

£ kil = £ j £ ak = E j { £ a k - E 

i=l j= 1 k=kj - i + l j = 1 k=l k= 1 

< £ a k + E j {<* - ( a - 

oo oo 

a +E^+E^°° 

3 = i i=i 

Lemma 3.3 Let X and Y be real separable Hilbert spaces with inner products 
denoted by < •,• >x and < •,• >y respectively. Then every $ e HS(X,Y ) can be 
written in factored form as $ = $ 1 $ 2 with e C(Y) compact and $ 2 e HS(X , Y). 

Proof Let be an orthonormal basis for X and let be an orthonormal 

basis for Y. If $ e HS(X , Y) then it has a representation in the form of an infinite 
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I 


matrix $ *-* [ <fij ] where ip tj =< y,, $Xj >y , and Y Y < 00 • Now for ' l = l> 2 v 

1=1 J = 1 

OO 

set a, = Y Vij- The sequence {a,}°^j is absolutely summable, so applying Lemma 
i=i 

3.2 we obtain sequences and {c,}?^ for which a, = 6,c t -, i = 1,2,..., lim b t = 

2 — ► OO 

OO OO 00 

0, and X^ |cj | = X) c i < °°- Define <2>* e £(Y) by § l y = X] \/&7 < y,y* >y y*, y e Y, 


2=1 


2=1 


2=1 


and <3> 2 e £(X,Y) by $ 2 x = X) XI Tyfr < >x y *’ x 6 X - Then $ 2 $ 2 = $, 

2=1 j = 1 v 1 

OO OO / \ 2 00 00 00 00 

and, since lim y/bl = 0 and Y Y \VFJ = X) E <Pij = E ff = E < °°> 

t =i j=i v v ‘ ' »=i j = 1 i=i i=i 

it follows that is compact and <F 2 e HS(X,Y). 


Lemma 3.4 


(a) For every lim |P n $P n — = 0. 

n — ►oo 

(b) For every V, lim ||P n $P n — <F|| V = 0. 

71 — ^OO 

Proof (a). For $ e H we have 

|P n $P n - *\ n < | P„$P„ - P„$|* H + |P n $ - *\ n 
< |*P„ - 4|„ + |P„* - #|„ 

= |(*P„)* - **|„ + |P„<f - »|* 

= |P„«* - «*|„ + |P„* - *|„ 

where in the above estimate we have used the fact that P n e C(H ) with |P n | < 1 
and that |P n 'I , |ft < |P„| I’Ll?* < |^|-« for every ’F e 'H. If we apply Lemma 3.3 with 
X = Y = H to <F, <F* e H = HS(H, H ), then we obtain $ = S 1 # 2 , <F* = ($*) 1 ($*) 2 
with <F 1 ,(<F*) 1 e £(P) compact and ($*) 2 e HS(H,H). It follows that 

I ?»•** - *'l„ = |(P» - /)(0*)'(4 *) 2 | h < |(P„ - /)($')■ 1 1(*-) 2 |„ , 

and 

|P»*-*l«<|(P»-0* , ll* 2 l«. 

which together with assumnption (3.1) and the fact that <F X and (<F*) 1 are compact 
yield the desired result. 

(b). For <F e V = HS(V*,H) D HS(H, V) we have 

|Pn^Pn — $\hS(V*,H) — — H S(V m ,H) + l-^n^ — $| HS(V m ,H) 

— |P n | |<FP„ — $\hS(V,H) + — ^1 HS(V*,H) 

— l^-fn — $ltfS(V*,tf) + _ ^l//S(V*,H). 


I 
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Now $ e HS(V*,H ) implies that $* e HS(H,V*) and that ($P„)* = 

P*$* e HS(H , V*) where, recalling Lemma 3.1, P* = 7 P„ 7 _1 is the adjoint of the 
operator P„ considered as an element of C(V*). It follows that 

(3.4) |P n $Pn - $IhS(V*,H) < I Pn^* ~ ®*\hS(H,V) + ~ $1 HS(V,H) 

Lemma 3.3 with X = H and Y = V* together with Lemma 3.1 imply that the first 
term on the right hand side of the estimate (3.4) tends to zero asn-+ oo. Similarly, 
Lemma 3.3 with X = V* and Y = H together with assumption (3.1) imply that 
the second term tends to zero as n — ► oo as well. A similar argument to the one 
given above can be used to show that lim |P n $P n — ^\hs(H V) = 0 an d the l emma 

is proved. 

The primary result of this paper is given in the following theorem. 

Theorem 3.1 Let II e V be the unique solution to the equation (2.5) (equivalently 
(2.3)) and for each n = 1,2,... let II n e 'H n be the unique solution to the approxi- 
mating operator equation (3.3). Then lim ||II n — II|| V = 0. 

n — *00 

Proof From the triangle inequality we obtain 


l|n„ - n|| v < ||n B - p„np n || v + ||p n np n - n|| v . 

An application of Lemma 3.4(b) yields that the second term on the right and side 
of the above estimate tends to zero as n — ► oo. As for the first term, we recall (2.4) 
and consider the estimate 

w p„ - p n np n || 2 v < [^{n„ - p n np n },n n - p„np n ]« 

= [An n + s„(n n ) -ail- B{ n), n n - p„np n ]* 

+ [^4II — *4P ^IlPrt, II n — P nJLP[ 

+[P(n) - H(p„np n ),n„ - p n up n ] n 

+ [B(P n UP n ) - Bn(U n ), U n - PnTLPnU 

= [©n ©5 IIn P nil P n]ft 

+[A{n - p n np n }, n n - p„np n ]* 

+[B( n) - s(p n np n ), n n - p n np n ]* 

-[p n p(n n )Pn - p(p r»np n), n„ - P„UP n ]n 

= [P n ©P„ - 0 ,n„ - PnTLPnh 
+[A{IL - p„np„}, n„ - P n UP n ]n 
+[ 5 (n) - 5(p n np n ),n n - p n np„]* 

~[B(u n ) - p(p„np n ),n n - pjip n ) n 

< [p n ©p„ - ©, n n - p n np n ] w 
+{A{ n - p n np n }, n„ - p„iip„]* 
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+[£(11) - £(P„IIP n ),n n - P n UP n ) n 

where in the final estimate above we have applied assumption (2.1). Continuing we 
find that 

w ||n„ - p n np n ||y < II PnQPn - ©Hv ||n„ - p n np n || v 

+ l-'^l^y V*) 11^ PnllP ally ll^n P rcUP n||y 

+ ||B(n)- p(p„np n )|| v . ||n„ - p„np n || v , 

or, recalling the continuous embedding of H into V* , that 

||n„ - p„np n || v < k | p n 0P n - ©|„ + \A\ c( w) lin - Pnnp n || v 

+K \B(Tl) - B(P n UP n )\ n . 

Lemma 3.4 together with the continuity assumption on B imply 

lim ||n„ - P n IIP n || v = 0 

n — ► oo 

and the theorem is established. 


4. An Example 

In order to illustrate the application of our thoery we consider an operator alge- 
braic Riccati equation arising in the design of an optimal feedback control law for a 
one dimensional heat/diffusion equation. Let H,U = 1/2(0, 1), both endowed with 
the usual inner product, (ip, tp) = f* ip(ij)tp(r})djj, and consider the linear-quadratic 
optimal control problem of finding u e L 2 (0,oo; U) which minimizes the quadratic 
performance index 



(Qx(t, •), x(t, •)) + r(u(t, •), u(t, ■ ))dt 


subject to the linear dynamical system 

(4- 1 ) §f(<,»7) “ = bu(t,ri),t > 0,0 < 77 < 1, 

(4.2) x(t, 0) = 0, x(t, 1) = 0, t> 0 

(4.3) x(0,r}) = x°(ri), 0 < rj < 1, 

where Q is a self-adjoint, nonnegative, and Hilbert-Schmidt operator from L 2 (0, 1) 
into P 2 (0,l),a e 1/00(0,1), a(r 7) > a > 0, a.e. 77 e (0,1), 6 e R, and x° e L 2 (0,1). 
If we define C = {$ e P5'(I/ 2 (0, 1), L 2 (0, 1)), $ = $*,$ > 0}, then Q e C and 

(Q<p)(r)) = Jo q(r], C)v>(C)^C with Q £ ^((0, 1) X (0, 1)), q(rj, Q = q((, 77), q(r), Q > 0, 
a.e. (77, C) e (0,1) x (0,1). 

Define V — Pq( 0, 1) endowed with the standard inner product, < ip, xp >= 
f 0 D<p(r})Dtp(r))dr} and corresponding norm, || • ||. It follows that V is densely, con- 
tinuously, and compactly embedded in H, that V* = H~ 1 (0, 1), and that H is 
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densely and continuously embedded in V *. Define the operator A e C(V,V*) by 
(A<p,rp) = (aD<p,D%l>), for ip,rj> e V. It follows that ( Aip,<p ) > c*j| 9 ?|| 2 , ip e V, and 
that the restriction —A of the operator -A to the set Dom(A) = {<p e V : A<p e H} 
(= H 2 ( 0, 1) fl Hi (0, 1) when a is suffciently smooth) is densely defined in H, nega- 
tive, self-adjoint, and it is the infinitesimal generator of a uniformly exponentially 
stable analytic semigroup, {T(t) : t > 0} of bounded, self-adjoint linear opeators on 
H. The solution to the initial-boundary value problem (4.1)-(4.3) is given by 

x(t ) = T(t)x° + f T(t — s)bu(s)ds, t > 0 
Jo 

where for each t > 0 x(t ) = x(t, •) e H = Z^O, 1) and u(t) = u(t, •) e U = ^2(0, 1) 
for almost every t > 0. 

The solution to the optimal control problem (see [5]) is given in closed-loop linear 
state feedback form by 


u(t) = — (b/r)Ux(t), a.e. t > 0 

where II is the unique nonnegative self-adjoint solution to the algebraic Riccati 
equation 

(4.4) A*n + nA + (6 2 /r)n 2 = 

It is clear that the existence-uniqueness and approximation theories presented above 
apply with 0 = Q e C and #($) = (6 2 /r)$ 2 for $ e C. It follows that there exists a 
unique solution II e HS(L 2 ( 0 , 1), ^ 2 ( 0 , 1)) to the nonlinear operator equation (4.4) 
with n = n*,n > 0, and n e ff5(Z 2 (0,l),ff 0 1 (0,l))nir5(Zr 1 (0,l),Z2(0,l)). Re- 
calling that HS(L 2 (0, 1), Z 2 ( 0, 1)) is isometrically isomorphic to Z 2 ((0, 1) x (0, 1)), 
II e C implies that there exists 7r e Z 2 ((0,1) x (0,1)) with 7 t(t 7,C) = ^(C? 7 ?) an( i 
7r(r7, C) > 0 for almost every (??, () e (0, 1) x (0, 1) for which 

= -(b/r) f ?r(?7,C)x(<,C)dC, 

Jo 

for almost every rj e (0, 1) and t > 0. 

We introduce linear spline based approximation. For each n=2,3,... let H n = 
span where for j=l,2,...,n-l, <p 3 n denotes the jth standard linear spline func- 

tion defined on the interval [0,1] with respect to the uniform mesh {0, 1/n, 2/n, ..., 1}. 
More precisely, 

0 < rj < j~ 1 /n 

j- 1 < r) <j/n 
j/n <7? <j + 1/n 
j + 1/n < rj < 1, 


^n(»7) 




0 

nr] - j + 1 
j A l - nr] 
0 
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j=l,2,...,n-l. Let P n : H —> H n denote the orthogonal projection of H onto H n 
with respect to the usual inner product on H = £ 2 ( 0 , 1) and define Galerkin ap- 
proximations A n e C(H n ) to A in the usual way. That is let A n (p n = ifi n where 
for ip n e H n ,tp n is the unique element in H n which satisfies (Aip n ,x n ) = (i’mXn), 
Xn e H n . Set Q n = P n Q e C(H n ). 

Using well known approximation properties of linear interpolatory splines (see [7]) 
it is not difficult to argue that lim \\P n <p — v?|| = 0, v? e Hq( 0,1) and consequently 
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that assumption (3.1) is satisfied. It follows therefore, from the theory presented 
in section 3 above, that there exists a unique nonnegative self-adjoint operator 
n n e C(H n ) satisfying the algebraic Riccati equation in H n given by 

(4.5) A*n n + U n A + (b 2 /r)H 2 n = Q n . 

Moreover, we have the Hilbert-Schmidt norm convergence of II n P„ to II as n — > 00 . 
That is 

( 4 * 6 ) | n „ p „ — n | hs(h,h) = 

We in fact also obtain that lim |II n P n — y) — 0 and that 

— V* ,H) = 0- 

From a computational point of view, since the basis elements (f 3 n are not mu- 
tually orthonormal, simply replacing the operators in the finite dimensional alge- 
braic Riccati equation (4.5) with their corresponding matrix representations will 
not lead to the usual symmetric matrix Riccati equation for which a variety of 
computational solution techniques exist. Toward this end, for a linear opera- 
tor L n with domain and/or range in H n , we denote its matrix representation 
with respect to the basis {<fi 3 n } n = l for H n by Ln- Define 4> n : [0,1] — * |R n_1 by 
= (^»(*7)»¥>n(*/)»--- iVn~ 1 ( r l)) T and set M N = (<£„,$[) = 
fo i)dn- If then follows that An = M^ 1 (aD<f> n ,D(f>n) w ith 

A* n = M ^ 1 AJjMn, and Qn = M^ i {Q(j) n ^(f>^). If we set Qn — MnQn and 
U N = M N n N , then II yv is the unique nonnegative self-adjoint solution to the 
(n — 1) x (n — 1) matrix algebraic Riccati equation given by 

(4.7) A^Un + IL n A n + (b 2 /r)fl N M^fl N = Q N . 

The approximating optimal control laws take the form 

u n (t,v) = ~(b/r) [ 7r n (?7,CMt,C)dC 

JO 
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for almost every r) t (0, 1) and t > 0 where 

*„(!»,<) = M’>) T MH , n N M?MQ, 

(rj,Q e [0,1] x [0,1]. It follows from (4.6) that lim 7r n = n in I/2((0, 1) x (0,1)) - 

n — >-oo 

that is that 

(4.8) lim fg f 1 l7r n (rf,C)-7r(Tj,Ol 2 d(dr}. 

n— +oo 

To illustrate we take a(r/) = a > 0, a constant, and let Q e HS(L2(0, 1), £2(0, 1)) 
be the finite rank modal projection operator given by 

V 

Qv = <P e L 2 (0,1) 

k = 1 


where v < 00, if’k( T l) = V^sink-Kr], r) e [0,1], k = 1,2, ...,v, and 9* > 0, k — 
1,2, A somewhat tedious, but rather straight forward computation yields 

V 

(Qn)j-; = ^ 2 qkSkiSkj, i,j = 1,2 , ..., n — 1, 
k = 1 


where 


5S, - (sin k h£±ll - 2 sin k*± + sin , 

Kl (kn) 2 ( n n n J 


k = 1,2, . . . ,u,£ = 1,2, ...,n — 1. Setting a = .25, b = 1.0, r = .01, v = 3, 
and 91 = 92 = 93 = 1.0, and using Schur- vector decomposition of the associated 
Hamiltonian matrix (see [4]) to solve the matrix Riccati equation (4.7) for various 
values of n we obtained the kernels, 7r n , plotted in the figures below. That the 
convergence given in (4.8) above is achieved is immediately clear. 
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